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We present a recent result on the existence of the dynamics in the thermodynamic limit 
of a class of anharmonic quantum oscillator lattices, which was obtained using Lieb— 
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•/^ , 1. Introduction 

OO ' In. condensed matter physics, three common types of degrees of freedom are often 

OA I found together: atoms positioned in a lattice, spin magnetic moments, and itinerant 

0^ ' electrons. Many phenomena primarily involve only one of these and it has therefore 

, proved useful to study them with separate models. Here, we are concerned with the 

spatial degrees of freedom of atoms in a lattice. The fact that the atoms oscillate 
about their equilibrium positions is essential for many important phenomena in 
condensed matter physics. The question we address is whether one can define the 
time evolution, consistent with the Schrodinger equation, for an infinite assembly 
^ . of quantum oscillators such as a crystal lattice. 

It is natural to start by considering the standard harmonic interaction, which 
should be an accurate description when the displacements from the equilibrium 
positions are small. It is well known, however, that some basic phenomena require 
that we consider anharmonic perturbations. Our approach applies to multi-body 
interactions that fall off sufficiently fast but for clarity and space limitations we 
will restrict ourselves to local and nearest neighbor anharmonicities here and refer 
to our forthcoming paper for the general case.^^ We will also limit the discussion 
to oscillators organized in a translation-invariant fashion on a lattice. To keep the 
notation simple we will work with one-dimensional oscillators at each lattice site, 
but this is not essential. With some straightforward modifications of the analysis 
our main result can be generalized to non-translation invariant models defined on 
a graph satisfying a few natural assumptions. 
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2. Harmonic lattice systems 

Let A denote a finite subset of tiie i^-dimensional hypercubic lattice iy and define 
Hamiltonians li\ as self-adjoint operators on Ha = (^) by the following 

expression 

where o), A, > 0, and px and qx are the canonical monicntum and position operators 
for the oscillator at a; € A. The Heiscnbcrg dynamics, {rj'^lfgR, is defined by 

The algebra of observables is a tensor product 

B{nK) = (^B{nx) = AK 

xeA 

such that for X c A, we have Ax C Ak, by identifying A e Ax with A^lf^\x € ^a- 
Therefore, for A e Ax, '''{^{A) is well-defined for all A containing X. 

Our main concern is the existence of a limiting dynamics: do we have Tt such 
that 

T^(A) n{A), as A T 

in a suitable sense? For an anharmonic lattice of classical oscillators a positive 
answer to this question was obtained by Lanford, Lebowitz and Lieb.^ 

It is well-known that in the case of quantum spin systems (i.e., bounded lo- 
cal Hamiltonians) one can use Lieb-Robinson bounds to establish the existence of 
the thermodynamic limit. ^ The essential observation is as follows. Let A„ be an 
increasing exhausting sequence of finite volumes with Hamiltonians of the form 

XCA„ 

where '^{X) = $(X)* G Ax- Then, for n > m, one easily derives the bound: 

Therefore, if we can show that the commutators in the integrand have sufficiently 
small norms, it will follow that the finite- volume dynamics form a Cauchy sequence. 
Estimates for such commutators were first derived by Lieb and Robinson.^ For 
A e Ax and B e Ay, they proved a bound of the form 

where C, a, and v arc positive constants and d{X, Y) denotes the distance between 
X and Y. Estimates of this type are now commonly referred to as Lieb-Robinson 
bounds. For anharmonic lattice systems Lieb-Robinson bounds were recently 
proved in Ref. 10, and this work builds on the results obtained there. 
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For the oscillator lattices with A 0, the approach suggested by Eq. [l]does not 
work due to the unboundedness of {qx — lyY ■ But we note that for the harmonic 
lattice, can be calculated explicitly and can also be defined for the infinite 
lattice. One then sees, however, that t{^{A) cannot converge in norm. The idea is to 
consider the infinite anharmonic system as a limit of finite- volume perturbations of 
the infinite harmonic system. We mention that, with a different approach. Amour, 
Levy-Bruhl, and Nourrigat have recently obtained convergence results for certain 
models by introducing suitable Sobolev-like norms for the observables.^ 

Rigorous perturbation theory of infinite systems is available if the unperturbed 
infinite system dynamics has a suitable continuity property.^ In our case this will 
be continuity for the weak operator topology, which follows from known properties 
and the exact solution of the harmonic lattice. 

Up to a redefinition of the parameters, the harmonic lattice model Hamiltonian 
on a finite subset A C Z'^, is 

acting on TL\. The creation and annihilation operators are defined by 

ttx = ^ {qx + iPx) and a* = {qx - ipx) , 

which satisfy the Canonical Commutation Relarions (CCR): 

[ax,ay] = [a*x, a*] = and [a^., a*] = Sx,y for aU x, y e A^ 

H\ is a quadratic expression in ax,a*, which can be diagonalized by a Bogoliubov 
transformation.® This is usually done in Fourier space resulting in: 

Ha^Y. 7(fc)(2fo^6fc + l) 

where 7(fc) = (w^ + 4AX]J=i sin^(fcj/2))^/^. "^^^ K ^""^ creation operators for the 
eigenmodes of the system, which also satisfy the CCR, and the ground state of 
H\ is the corresponding vacuum. The time evolution is simply given by Tt{h*^) = 
g-2i7(fc)t^* ^ It is useful for us to express the Bogoliubov transformations in real 
space: for each / : A^ — > C, set 

«(/) = E /(^«- ^ E /(^x- 

and similarly for the Vs. Then, there are real-linear operators U and V such that 

a(/) = b{U*f) - b*{V*f), a*{f) - b*{U*f) ~ b*{V*f). 
Using these relations, one finds a one-parameter group, Tt, such that 



rMf)) = n{b{U*f) - b*{V*f)) = a{TJ), n{a*{f)) = a*(TJ). 
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Since the canonical operators (p's, q's , a's, a*'s,...) are all unbounded, we will derive 
Lieb-Robinson bounds for the Weyl operators instead: 



— {a{f) + a*{f)) 



Wif) = exp 

This can be done very explicitly by using the relations 

Tt{W{f)) = W{Ttf), and W{f)W{g) = e'^^^'^) W(/ + g) 
where (T{f,g) = lm{f,g). Observe that 

[rtmf)),W{g)] = {W{TJ) - W{g)WiTJ)W{-g)}W{g) 
= |l _ e-2-(^*/'5)} W(Tt/)W(5) . 

Since the Weyl operators are unitary, we therefore have ||[Tf(VF(/)), VF((;)]|| < 
2\a{Ttf,g)\. So, all we have to do is to estimate a{Ttf,g). Ttf is explicitly given by 

^t/ = / * h[^t + / * ^2^t' where the functions hi^t and /i2,t are given by 



keAi 



feeA* 



fceA? 



By analyzing these functions one then shows that for every a > 0, there exist Ca 
and Va such that 



W{Ttf,g)\<Cae-^^'^J2\f{x)\\g{y)\e-^'^^' 



v) 



Theorem 2.1 (Ref. 10). Let X,u) > 0. Then, for all f,g with supp/ C X and 
suppg C Y, 

[rt{W{f)),W{g)] II < Cj2\m\ l5(2/)|e-'^'^"'^'-''l*l^ 
x,y 

with V = Qy/co'^ + AvX. 



3. Anharmonic lattice systems 

In Ref. 10 we proved Lieb-Robinson bounds for finite systems, with estimates uni- 
form in the volume. With the same approach one can also prove bounds for harmonic 
infinite systems with an anharmonic perturbation in finite volume, i.e., a model with 
formal Hamiltonian of the form: 

Typical $ are i>({a;}) — V{qx), ^{{x,y}) = W{qx — Qy), etc. but more general 
perturbations can be considered. 
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For the harmonic system, we can directly see how the exact expressions extend 
to the infinite system. For / : C in a suitable function space (£^{'Zi'^), or 

^^(Z'')), it is straightforward to define 

Tt{Wif))^W{TJ) 

We get a Hilbert space representation by representing the Weyl operators on the 
Fock space generated by the b*{f) operators acting on the vacuum. On this space 
we have well-defined Hamiltonian such that 

Tt{W{f))^e'*''W{f)e-''" 

Tt is the dynamics for the formal Hamiltonian 

and the Lieb-Robinson bounds continue to hold. 

For simphcity, consider the perturbation of the form Pa = J2xeA^(lx)- Then, 
the perturbed dynamics, formally corresponding to H + Pa, and can be defined 
mathematically by the Dyson series: 



'A),[---[rtdPA),rt{W{m]] 



n=l J0<ti<t2-<t 

We have the following Lieb-Robinson bounds for t^. 

Theorem 3.1 (Ref. 10). Let A > 0,w > 0, and V such that \\PV{k)\\i < oo. 
Then, for all f,g (z t^{'L^), we have 

I [T^{W{f)),W{g)] II <C^|/(x)||g(y)|e-2('^(-^^)-''l*l) 

x,y 

with 



V = 6Vw2 + 4i/A + c||fc2T/(fc)||i. 



To show convergence, we estimate ||tj "{W{f)) — (T/F(/))||, for Am C A„ by 
considering t^" as a perturbation of t^"^ . This gives 



Therefore 



'{W{f)) 



26A„\A„ 



'\t\-s 



{W{f)) 



By writing V{qx) = J V{p)W{pSx)dp , W{p6x) — e^P'^" , we can then use Theorem 



EH for 



to obtain the convergence. The convergence is uniform on intervals 



[—to, to]- We also immediately get continuity in t of the limiting dynamics by an e/3 
argument. 



September 15, 2009 17:57 



WSPC - Proceedings Trim Size: 9.75In x 6.5in icmp'prague 



6 

Theorem 3.2 (Ref. 11). Assume u; > 0, \\kV{k)\\i < oo,\\k'^V{k)\\i < oo. For 
all f and all t G R, the limit 

converges in the operator norm topology and the resulting the dynamics is continuous 
in t in the weak operator topology. 
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